function [Location,DistMatrix,Ncities,x,val] = pr76init(InitOps); 
Location = [3600 2300;
3100 3300;
4700 5750;
5400 5750;
5608 7103;
4493 7102;
3600 6950;
3100 7250;
4700 8450;
5400 8450;
5610 10053;
4492 10052;
3600 10800;
3100 10950;
4700 11650;
5400 11650;
6650 10800;
7300 10950;
7300 7250;
6650 6950;
7300 3300;
6650 2300;
5400 1600;
8350 2300;
7850 3300;
9450 5750;
10150 5750;
10358 7103;
9243 7102;
8350 6950;
7850 7250;
9450 8450;
10150 8450;
10360 10053;
9242 10052;
8350 10800;
7850 10950;
9450 11650;
10150 11650;
11400 10800;
12050 10950;
12050 7250;
11400 6950;
12050 3300;
11400 2300;
10150 1600;
13100 2300;
12600 3300;
14200 5750;
14900 5750;
15108 7103;
13993 7102;
13100 6950;
12600 7250;
14200 8450;
14900 8450;
15110 10053;
13992 10052;
13100 10800;
12600 10950;
14200 11650;
14900 11650;
16150 10800;
16800 10950;
16800 7250;
16150 6950;
16800 3300;
16150 2300;
14900 1600;
19800 800;
19800 10000;
19800 11900;
19800 12200;
200 12200;
200 1100;
200 800];
Ncities= size ( Location, 1) ;
DistMatrix= dists ( Location, Location) ;
x=[1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76];
x=x';
val=0;
for i=1:Ncities-1
    val=val+DistMatrix(x(i),x(i+1));
end
val=val+DistMatrix(x(1),x(Ncities));


